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ABSTRACT 

We present the discovery of five quasars at z ~ 6 selected from 260 deg 2 of the Sloan Digital 
Sky Survey (SDSS) southern survey, a deep imaging survey obtained by repeatedly scanning a stripe 
along the Celestial Equator. The five quasars with 20 < zab < 21 are 1-2 magnitudes fainter than 
the luminous z ~ 6 quasars discovered in the SDSS main survey. One of them was independently 
discovered by the UKIRT Infrared Deep Sky Survey. These quasars, combined with another z ~ 6 
quasar known in this region, make a complete flux- limited quasar sample at zab < 21. The sample 
spans the redshift range 5.85 < z < 6.12 and the luminosity range —26.5 < M1450 < —25.4 (Ho = 70 
km s _1 Mpc -1 , fl m = 0.3, and Q\ = 0.7). We use the 1/V a method to determine that the comoving 
quasar spatial density at (z) = 6.0 and (A/1450} = -25.8 is (5.0 ± 2.1) x 10~ 9 Mpc~ 3 mag -1 . We 

model the bright-end quasar luminosity function (QLF) at z ~ 6 as a power law $(£1450) oc L^ 450 . 
The slope [3 calculated from a combination of our sample and the luminous SDSS quasar sample is 
—3.1 ± 0.4, significantly steeper than the slope of the QLF at z ~ 4. Based on the derived QLF, 
we find that the quasar/AGN population cannot provide enough photons to ionize the intergalactic 
medium (IGM) at z ~ 6 unless the IGM is very homogeneous and the luminosity (L| 450 ) at which 
the QLF power law breaks is very low. 

Subject headings: galaxies: active — quasars: emission lines — quasars: general 



1. INTRODUCTION 

High-redshift quasars are among the most luminous 
objects known and provide direct probes of the distant 
universe when the first generation of galaxies and quasars 
formed. In recent years, over twenty z ~ 6 quasars with 
zab < 20 have been discovered (e.g. Fan et al. 2000, 
2001a, 2003, 2004, 2006a; Goto 2006). These luminous 
quasars are essential for understanding the accretion his- 
tory of black holes (BHs), galaxy formation, and chemical 
evolution at very early epochs. They harbor supermas- 
sive BHs with masses higher than 10 9 M@ and emit near 
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the Eddington limit (e.g. Barth et al. 2003; Vestergaard 
2004; Jiang et al. 2006a; Kurk et al. 2007), revealing 
the rapid growth of central BHs at high redshift. Their 
emission lines show solar or supersolar metallicity in the 
broad line regions, indicating that there was vigorous star 
formation and element enrichment in the first gigayear of 
cosmic time (e.g. Barth et al. 2003; Maiolino et al. 2003; 
Jiang et al. 2007; Kurk et al. 2007). Their absorption 
spectra show that the intergalactic medium (IGM) at 
z ~ 6 is close to the reionization epoch (e.g. Becker et al. 
2001; Djorgovski et al. 2001; Fan et al. 2006b,c). 

The majority of the currently known z ~ 6 quasars 
were discovered from ~ 8000 deg 2 of imaging data of 
the Sloan Digital Sky Survey (SDSS; York et al. 2000). 
They were selected as z-dropout objects using optical col- 
ors. Several other high-redshift quasars were discovered 
based on their infrared or radio emission. For example, 
Cool et al. (2006) discovered one quasar at z = 5.85 in 
the NOAO Deep Wide-Field Survey (NDWFS; Jannuzi & 
Dey 1999) Bootes Field using the AGN and Galaxy Evo- 
lution Survey (AGES) spectroscopic observations. The 
quasar was selected from a Spitzer mid-infrared quasar 
sample and has a zab magnitude of 20.68 and an optical 
luminosity of M B = -26.52. By matching the FLAMIN- 
GOS Extragalactic Survey IR survey (Elston et al. 2006) 
data to the Faint Images of the Radio Sky at Twenty- 
cm (FIRST; Becker et al. 1995) data, McGreer et al. 
(2006) discovered a radio-loud quasar at z = 6.12 in 4 
deg 2 of the NDWFS region. This quasar is a broad ab- 
sorption line (BAL) quasar with an optical luminosity of 
Mb = —26.9, comparable to the luminous SDSS quasars 
at z ~ 6. 

Despite the high-redshift quasar surveys mentioned 
above, very little is known about faint quasars (zab > 
20) at z ~ 6. The SDSS main survey only probes the 
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most luminous quasars, and with a density of 1/470 deg 2 
(Fan et al. 2006a). The Cool et al. (2006) quasar at 
z = 5.85 was zab > 20, but the sample contains a single 
object and is selected from an area of less than 10 deg 2 . 
Mahabal et al. (2005) found a very faint quasar with 
zab = 23.0 at z = 5.70 in a 2.5 deg 2 field around the lu- 
minous quasar SDSS J114816.64+525150.3 11 at z = 6.42. 
Willott et al. (2005) imaged a 3.83 deg 2 region down to 
Zab = 23.35 in the first results of the Canada-France 
High-redshift Quasar Survey (CFHQS) and did not find 
any quasars at z > 5.7. In these surveys both the quasar 
samples and the survey areas are very small, thus they 
do not provide a good statistical study of high-redshift 
quasars at zab > 20. Recently, Willott et al. (2007) dis- 
covered four quasars at z > 6 from about 400 deg 2 of 
the CFHQS, including the most distant known quasar 
at z = 6.43. Three of these quasars have zab magni- 
tudes fainter than 21. Since their follow-up observations 
are not yet complete, they did not determine the spatial 
density of these quasars. 

Finding faint quasars at z ~ 6 is important for study- 
ing the evolution of the quasar population and quasars' 
impact on their environments. Fan et al. (2004) ob- 
tained the bright-end quasar luminosity function (QLF) 
at z ~ 6, but the slope, — 3.2±0.7, was very uncertain due 
to the small luminosity range of the sample. Richards et 
al. (2004) put a broad constraint on the bright-end slope 
of > —4.63 (3<T ) from the absence of lenses in four quasars 
at z ~ 6. Shankar & Mathur (2007) have considered the 
implications of all existing z ~ 6 quasar observations, 
including deep X-ray surveys, for the faint end of the 
high-redshift QLF. Based predominantly on the X-ray 
surveys, they argue that there is a flattening of the QLF 
at Mi4 50 > -24.67. 

With the discovery of faint high-redshift quasars, the 
QLF can be well determined. The QLF at z ~ 6 is im- 
portant to understand BH growth at early epochs (e.g. 
Volonteri & Rees 2006; Wyithe & Padmanabhan 2006). 
While bright quasars at high redshift have central BH 
masses between 10 9 and 10 10 M©, fainter quasars with 
zab > 20 are expected to harbor BHs with masses of a 
few times 10 8 M & or below (e.g. Kurk et al. 2007), which 
may be associated with galaxies of lower masses. The 
QLF also enables us to determine the quasar contribu- 
tion to the UV background at z ~ 6. Detection of com- 
plete Gunn-Peterson troughs (Gunn & Peterson 1965) 
among the highest-rcdshift quasars indicates a rapid in- 
crease of the IGM neutral fraction at z ~ 6, and suggests 
that we have reached the end of the reionization epoch 
(e.g. Becker et al. 2001; Djorgovski et al. 2001; Fan et 
al. 2006c). It is unclear what individual contributions of 
galaxies and quasars to the reionization are. Although 
there is evidence showing that quasars are probably not 
the main contributor to reionization (e.g. Salvaterra et 
al. 2007; Srbinovsky & Wyithe 2007; Shankar & Mathur 
2007), a proper determination of the QLF at z ~ 6 is 
needed to constrain the the quasar contribution. 

In this paper we present the discovery of five z ~ 6 
quasars with 20 < zab < 21 selected from 260 deg 2 of the 
SDSS southern survey, a deep imaging survey obtained 

11 The naming convention for SDSS sources is SDSS 
JHHMMSS.SSiDDMMSS.S, and the positions are expressed in 
J2000.0 coordinates. We use SDSS JHHMMiDDMM for brevity. 



by repeatedly scanning a 300 deg 2 area in the Fall Celes- 
tial Equatorial Stripe ( Adclman-McCarthy et al. 2007a) . 
One of the five quasars, SDSS J020332.39+001229.3 
(hereafter SDSS J0203+0012), was independently discov- 
ered by matching the UKIRT Infrared Deep Sky Survey 
(UKIDSS; Warren et al. 2007) data to the SDSS data 
(Venemans et al. 2007). These five quasars, together 
with another quasar, SDSS J000552. 34-000655.8 (here- 
after SDSS J0005-0006) previously discovered in this re- 
gion (Fan et al. 2004), form a well-defined low- luminosity 
quasar sample at high redshift. We use this sample and 
the luminous SDSS quasar sample to measure the QLF 
and constrain the quasar contribution to the reionization 
of the universe at z ~ 6. 

The structure of the paper is as follows. In § 2 we 
introduce the quasar selection criteria and photometric 
and spectroscopic observations of quasar candidates. In 
§ 3 we describe the properties of the five new quasars. 
We derive the QLF at z ~ 6 in § 4, and discuss the 
contribution of quasars to the ionizing background in § 
5. We give a brief summary in § 6. Throughout the paper 
we use a A-dominatcd flat cosmology with Ho = 70 km 
s" 1 Mpc -1 , fi m = 0.3, and fl A = 0.7 (Spergel et al. 
2007). 

2. CANDIDATE SELECTION AND OBSERVATION 
2.1. SDSS Deep Imaging Data 

The SDSS is an imaging and spectroscopic survey of 
the sky (York et al. 2000) using a dedicated wide-field 
2.5 m telescope (Gunn et al. 2006) at Apache Point Ob- 
servatory. Imaging is carried out in drift-scan mode us- 
ing a 142 mega-pixel camera (Gunn et al. 1998) which 
gathers data in five broad bands, ugriz, spanning the 
range from 3000 to 10,000 A (Fukugita et al. 1996), on 
moonless photometric (Hogg et al. 2001) nights of good 
seeing. The effective exposure time is 54 seconds. The 
images are processed using specialized software (Lupton 
et al. 2001), and are photometrically (Tucker et al. 2006; 
Ivezic et al. 2004) and astrometrically (Pier et al. 2003) 
calibrated using observations of a set of primary stan- 
dard stars (Smith et al. 2002) on a neighboring 20-inch 
telescope. All magnitudes are roughly on an AB system 
(Oke & Gunn 1983), and use the asinh scale described 
by Lupton ct al. (1999). 

A primary goal of the SDSS imaging survey is to scan 
8500 deg 2 of the north Galactic cap (hereafter referred 
to as the SDSS main survey). In addition to the main 
survey, SDSS also conducts a deep survey by repeat- 
edly imaging a 300 deg 2 area on the Celestial Equa- 
tor in the south Galactic cap in the Fall (hereafter re- 
ferred to as the SDSS deep survey; Adelman-McCarthy 
et al. 2007a). This deep stripe (also called Stripe 82) 
spans 20 h < RA < 4 h and -1.25° < Dec < 1.25°. 
The multi-epoch images, when coadded, allow the se- 
lection of much fainter quasars than the SDSS main sur- 
vey. Jiang et al. (2006b) have used the deep data to find 
low-redshift faint quasars selected from the SDSS coad- 
ded catalog, i.e., each run goes through the photometric 
pipeline PHOTO separately, and the resulting catalogs 
are coadded. The SDSS deep survey will eventually reach 
a depth of %ab = 24.0 ~ 24.5 (5 a detection for point 
sources) with more than 50 epochs of data (Adclman- 
McCarthy et al. 2007b). The area and the depth will 
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then be comparable to the CFHQS, which has covered 
~ 400 deg 2 to a limit of iab ~ 24.5 (Willott et al. 2007). 

At the time the coadded images used in this paper 
were made, 2005, a given area of sky on Stripe 82 had 
been scanned 5-18 times under standard SDSS imaging 
conditions, and all available data were included in this 
version of the coadds. The construction of the coadds is 
summarized as follows. First, each input image was cal- 
ibrated to a standard SDSS zeropoint using the median 
of the SDSS photometric solutions for the runs covering 
the area. For each image we made an inverse variance 
image to serve as a weight map, and within the map 
we assigned near-zero values to pixels with the INTERP 
bit set in the fpMask files of PHOTO, which indicates 
that a bad pixel, bad column, or cosmic ray has been 
interpolated over. After the images were sky subtracted 
using the PHOTO sky estimate, each image was mapped 
onto a uniform rectangular output astrometric grid using 
a modified version of the registration software SWARP 
(Bcrtin ct al. 2002). The main modification was to in- 
corporate the known SDSS camera distortions into the 
astromctry. The weight maps were subjected to the same 
mapping. The mapping used the LANCZOS3 kernel for 
the images and weights, while the INTERP bit images 
were mapped using NEAREST before being used to set 
near-zero values in the weight map images. The mapped 
images were then coadded using a flux scaled inverse- 
variance weighted average. The resulting images had 
known photometric and astrometric calibrations. For 
precision photometry, the PSFs as a function of posi- 
tion on the images were required. We used an algorithm 
(implemented in PHOTO) that coadded the PSFs as a 
function of position known from each input image using 
the same weights as the images were averaged. Once the 
photometry, astrometry, and PSF were known, we pro- 
ceeded to run a version of FRAMES (the main portion 
of PHOTO) slightly modified to take effective gains and 
sky noises into account when calculating error estimates. 
The uncalibrated output of PHOTO was adjusted and 
placed on the SDSS AB system using a slightly modifed 
version of the SDSS target selection software TARGET. 
The details of the construction of the coadds will be left 
for a forthcoming paper. 

The resulting data were tested in three ways. First, 
object by object comparisons were made against single- 
run SDSS data. Second, object by object comparisons 
were made with the SDSS coadded catalogs produced 
by suitably averaging all the catalog information from 
the individual runs covering Stripe 82 (courtesy of R. 
Scranton and D. Johnston). Third, statistical internal 
measurements, such as stellar color-color diagrams, were 
made. These tests show that the images were properly 
coadded and the depth of the coadds is close to what 
was expected from standard error propagation. Figure 
1 compares the photometric errors of PSF magnitudes 
for point sources in the SDSS main survey data and the 
SDSS coadded data. The error estimates were produced 
by the SDSS photometric pipeline PHOTO. From Figure 
1 the photometric errors in the coadded data are signifi- 
cantly smaller than those in the single-run data. 

In this paper we used the data in the range 310° < 
RA < 60°, as there were significantly fewer than 10 runs 
covering the range 300° < RA < 310°. The data also 
contains some "holes" in which the coadded images were 
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Fig. 1. — Photometric errors of PSF magnitudes for point sources 
as a function of %ab and zab magnitudes. The error estimates were 
produced by the SDSS photometric pipeline PHOTO. The errors 
in the SDSS coadded data are significantly smaller than those in 
the SDSS main survey data. They are close to those expected from 
1/2 

""singlc/^epochi where <T sing i e is the error in SDSS single-run data 
and iVepoch is the number of the runs used in the coadds. The flux 
limits of the z ~ 6 quasar surveys used in the SDSS main survey 
(e.g. Fan et al. 2001a) and in our work are indicated as the two 
vertical dotted lines, at which o(zab) ~ 0.1. 

not available. The effective area for this work is 260 deg 2 . 
The median seeing as measured in the riz bands was 
1.2" ± 0.05", where the error is the standard deviation 
of the seeing measured by PHOTO across the coadded 
images on Stripe 82. 

2.2. Quasar Selection Procedure 

Because of the rarity of high-redshift quasars and over- 
whelming number of contaminants, our selection proce- 
dure of z > 5.7 faint quasars from the multi-epoch SDSS 
imaging data contains the following separate steps (see 
also Fan et al. 2001a, 2003). 

1. Select i-dropout sources from the SDSS deep 
stripe. Objects with iab — zab > 2.2 and 
zab < 21 that were not detected in the ugr bands 
were selected as i-dropout objects. We rejected 
sources with one or more of the following SDSS 
processing flags: BRIGHT, EDGE, BLENDED, 
SATUR, MAYBE.CR, and MAYBE.EGHOST 
(see Stoughton ct al. 2002). At z > 5.7, the Lya 
emission line begins to move out of the SDSS i fil- 
ter, so a simple cut of iab — zab > 2.2 is used 
to separate high-redshift quasars (and cool brown 
dwarfs) from the majority of stellar objects (e.g. 
Fan et al. 2001a). At zab — 21, the photomet- 
ric errors of the coadded data reach o~(zab) ~ 0.1 
as shown in Figure 1, so the zab < 21 criterion 
guarantees a high quasar selection efficiency due to 
small photometric errors. 

2. Remove false i-dropout objects. All i-dropout ob- 
jects were visually inspected, and false detections 
were deleted from the list of candidates. The ma- 
jority of the contaminants are cosmic rays. Al- 
though the SDSS photometric pipeline effectively 
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Fig. 2. — The i\B ~ Z AB vs - z AB ~ J (and H) color-color diagrams. The open circles represent known L/T dwarfs from Golimowski et al. 
(2004), Knapp ct al. (2004), and Chiu ct al. (2006). The crosses represent simulated quasars at 5.8 < z < 6.6. The median track of quasar 
colors as a function of redshift is shown as the filled circles. We use the dashed lines to separate high-rcdshift quasar candidates from L/T 
dwarfs. Note that in practice there is no unambiguous separation between the dwarf locus and the quasar locus due to photometric errors. 



rejects cosmic rays, the leakage of a tiny fraction of 
cosmic rays will contribute a large contamination 
to our sample. Cosmic rays were recognized by 
comparing the individual multi-epoch images mak- 
ing up the coadds. Brown dwarfs with high proper 
motions can be removed in a similar way, since the 
multi-epoch images were taken over a period of five 
years. In the selection of luminous z ~ 6 quasars 
in the SDSS main survey, Fan et al. (2001a) used 
an additional step, z-band photometry of i-dropout 
objects, to eliminate cosmic rays and improve the 
photometry of potential candidates. In this work 
we did not use this step, as the photometry of the 
coadds is robust. About 60 i-dropout objects re- 
mained in this step. 

3. Near- infrared (NIR) photometry of i-dropout ob- 
jects. We then carried out NIR (J or H band) 
photometry of j-dropout objects selected from the 
previous step. The details of the NIR observa- 
tions are described in § 2.3. Using the iab — zab 
vs. zab — J (or H) color-color diagrams (Fig- 
ure 2), high- redshift quasar candidates were sep- 
arated from brown dwarfs (L and T dwarfs), which 
have more than ten times higher surface density. 
The open circles in Figure 2 represent known L/T 
dwarfs from Golimowski et al. (2004) , Knapp et al. 
(2004), and Chiu ct al. (2006). The crosses repre- 
sent simulated quasars at 5.8 < z < 6.6. Although 
there is no clear separation between the dwarf locus 
and the quasar locus due to photometric errors, we 
selected quasars with the following criteria, 

iAB-ZAB > 2.2 && zab- J < 0.5(i A B-ZAB) + 0-5, (1) 

or %ab - zab > 2.2 && Z A B - H < 0.5(iab - zab) + 1-0. 

(2) 

A total of 14 objects satisfied the criteria. 

4. Follow-up spectroscopy of quasar candidates. The 
final step is to carry out optical spectroscopic ob- 



servations of quasar candidates to identify high- 
redshift quasars. The details of the spectroscopic 
observations are described in § 2.3. 

2.3. NIR Photometry and Optical Spectroscopic 
Observations 

In the first two steps we described above, we selected 
about 60 i-dropout objects with iAB — zab > 2.2 and 
zab < 21 from 260 deg 2 of the SDSS coadded imaging 
data. After the selection of i-dropouts, we carried out 
J or _ff-band 12 photometry of these i-dropouts using the 
SAO Widefield InfraRed Camera (SWIRC) on the MMT 
in November 2005 and October 2006 to separate quasar 
candidates and cool dwarfs. We used a 5 x 5 dither pat- 
tern to obtain good sky subtraction and to remove cosmic 
rays. The exposure time at each dither position was 30 
seconds. The total exposure time for each target was cal- 
culated to achieve an uncertainty of uj (or an) ~ 0.08. 
The typical exposure time on individual targets is 10 
minutes. The SWIRC data were reduced using standard 
IRAF 13 routines. Five bright objects from the Two Mi- 
cron All Sky Survey (2MASS; Skrutskie et al. 2006) in 
the field of each target were used to apply the aperture 
correction and absolute flux calibration. 

After the NIR photometry of the i-dropouts, we se- 
lected quasar candidates that satisfied the criteria of 
Equation 1 or 2. Optical spectroscopy of the candidates 
was carried out using the Echelle Spectrograph and Im- 
ager (ESI; Sheinis et al. 2002) on the Keck-II in Jan- 
uary 2006 and the Low Dispersion Survey Spectrograph 
(LDSS-3) on Magellan-II in October 2006. The observa- 
tions on ESI were performed in echellette mode, which 

12 The J-band photometry is more efficient for selecting z ~ 6 
quasars. Due to an instrument problem, however, the filter wheel 
was stuck and only the H filter was available on the night of Oc- 
tober 2006. 

13 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 
Optical and NIR Photometry 



Quasar (SDSS) 



Redshift a %ab (mag) zab (mag) J (mag) 



H (mag) 



J000552.34-000655.8 b 
J020332.39+001229.3 C 
J030331.40-001912.9 
J035349.72+010404.4 
J205406.49-000514.8 
J231546.57-002358.1 



5.850±0.003 
5.854±0.002 
6.070±0.001 
6.049±0.004 
6.062±0.004 
6.117±0.006 



23.40 ± 0.34 
23.72 ± 0.22 
23.92 ±0.23 
24.03 ± 0.30 
23.30 ±0.22 
24.90 ± 0.28 



20.54 ±0.10 

20.87 ±0.10 
20.85 ±0.07 
20.54 ± 0.08 
20.72 ± 0.09 

20.88 ± 0.08 



19.87 ±0.10 
19.05 ±0.08 



19.18 ±0.06 
19.94 ±0.08 



19.46 ±0.10 
18.55 ±0.06 



Note. — The %ab and zab magnitudes are AB magnitudes and the J and H magnitudes are 
Vega-based magnitudes. 

a The errors of the redshifts are the uncertainties obtained from our fitting process. 

b This quasar was discovered by Fan et al. (2004). The magnitudes were taken from Fan et al. 

(2004), and the redshift was determined from the Mgn emission line by Kurk et al. (2007). 

c This quasar was independently discovered by Venemans et al. (2007). 



provides excellent sensitivity from 4000 to 10,000 A. The 
observations on LDSS-3 were performed in longslit mode. 
LDSS-3 was designed to be very red sensitive. We used 
the VPH red grism with a ruling density of 660 lines/mm. 
The VPH red grism offers excellent throughput in the 
wavelength range from 6000 to 10,000 A. The exposure 
time for each target was 20 minutes, which is sufficient to 
identify z ~ 6 quasars in our sample. If a target was iden- 
tified as a quasar, several further exposures were taken 
to improve the spectral quality. The quasar data were 
reduced using standard routines. After bias subtraction, 
flat-fielding, and wavelength calibration were applied to 
the frames, one-dimensional spectra were extracted, and 
were flux calibrated using the spectra of spectroscopic 
standard stars. 

3. DISCOVERY OF FIVE NEW QUASARS AT 
Z ~ 6 



From the spectroscopic observations on Keck/ESI and 
Magcllan/LDSS-3 we discovered five z ~ 6 quasars in the 
SDSS deep stripe. Their z-band finding charts arc shown 
in Figure 3. Note that SDSS J0203±0012 was indepen- 
dently discovered by Venemans et al. (2007). Another 
z <~ 6 quasar in this area, SDSS J0005-0006 discovered 
by Fan et al. (2004) , was recovered by our selection crite- 
ria. These six quasars comprise a complete flux-limited 
sample at zab < 21. The optical and NIR properties of 
the quasars are given in Table 1. The iab and zab mag- 
nitudes of the newly discovered quasars are taken from 
the SDSS deep imaging data, and their J and H magni- 
tudes are obtained from our MMT/SWIRC observations. 
All the quasars have zab magnitudes between 20 and 21. 
The surface density of z ~ 6 quasars with zab < 20 is 
about 1/470 deg 2 (Fan et al. 2006a), so it is reasonable 
to find no quasars with zab < 20 in a 260 deg 2 area. 
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The optical spectra of the six quasars are shown 
in Figure 4. The spectrum of SDSS J0005-0006 was 
taken from Fan ct al. (2004). The spectra of SDSS 
J035349. 72+010404.4 (hereafter SDSS J0353+0104) and 
SDSS J231546.57-002358.1 (hereafter SDSS J2315-0023) 
were taken on Keck/ESI with a total exposure time of 
60 minutes on each source. The spectra of the other 
three quasars were obtained with Magellan/LDSS-3, and 
the total exposure time on each source was 100 minutes. 
Each spectrum shown in Figure 4 has been scaled to the 
corresponding zab magnitude given in Table 1, thereby 
placed on an absolute flux scale. 

We estimate the redshifts for the new quasars from 
either the Lya, Nv A1240 (hereafter Nv), or the Oi 
A1304 (hereafter Oi) emission line. For each quasar, 
we measure the line center of one strong emission line 
using a Gaussian profile to fit the top ~ 50% of the 
line. This provides a rough estimate of the redshift. 
Using this redshift we subtract the power-law contin- 
uum and decompose the blended Lya and Nv emission 
lines into individual components. The details are de- 
scribed in the next paragraph. Then wc determine the 
redshifts from individual emission lines. The redshift of 
SDSS J030331.40-001912.9 (hereafter SDSS J0303-0019) 
is measured from the N v emission line, which is well sep- 
arated from Lya due to the narrow line width. The mea- 
sured redshift 6.070±0.001 is consistent with the redshift 
6.069±0.002 determined from the weak O I emission line. 
SDSS J0353+0104 is a BAL quasar as seen from strong 
absorption features around Lya, so its redshift is mea- 
sured from the Oi emission line. The redshifts of the 
other three quasars are estimated from the Lya emission 
lines. They are usually biased because the blue side of 
Lya is affected by the Lya forest absorption (Schneider 
et al. 1991). The mean shift with respect to the systemic 
redshift at z > 3 is about 600 km s" 1 (Shen et al. 2007), 
corresponding to 6z ~ 0.015 at z ~ 6. We correct for this 
bias for the redshifts measured from Lya. The results are 
listed in Column 2 of Table 1. The errors in the table 
are the uncertainties obtained from our fitting process. 
For the redshifts measured from Lya, their real errors 
could be much larger due to the scatter in the relation 
between Lya redshifts and systemic redshifts (Shen et al. 
2007). In our sample four quasars have redshifts greater 
than 6. The most distant quasar, SDSS J2315-0023, is 
at z = 6.12. 

We measure the rest-frame equivalent width (EW) and 
full width at half maximum (FWHM) of Lya and N V for 
each quasar except the BAL quasar SDSS J0353+0104. 
To allow the analysis of the emission lines we first fit 
and subtract the continuum. The wavelength coverage 
of each spectrum is too short to fit the continuum slope, 
so we assume it is a power law with a slope a v = —0.5 
{fv ~ v a "), and normalize it to the spectrum at rest 
frame 1275-1295 A, a continuum window with little con- 
tribution from line emission. Lya and Nv are usually 
blended with each other, so we use three Gaussian pro- 
files to simultaneously fit the two lines, with the first 
two profiles representing broad and narrow components 
of Lya and the third representing Nv. Since the blue 
side of the Lya emission line is strongly absorbed by 
Lya forest absorption systems, we only fit the red side 
of the line and assume that the line is symmetric. We 



ignore the weak Sin A1262 emission line on the red side 
of N v. The measured EW and FWHM in units of A are 
shown in Table 2. We also give the FWHM of Lya in 
units of km s -1 . We emphasize that the EW and FWHM 
of Lya in the table have taken into account the absorbed 
emission by the Lya forest, while most previous studies 
did not take this absorption into account. 

The distributions of the Lya EW and FWHM are 
broad. The average Lya EW and FWHM measured 
from the low-redshift SDSS composite spectrum of Van- 
den Berk et al. (2001) are about 90 A and 20 A (~ 5000 
km s _1 ), respectively. The EW of Lya+N V from a sam- 
ple of quasars at 3.6 < z < 5.0 is 69±18 A (Fan et al. 
2001b), although this is affected by the Lya forest ab- 
sorption. We analyzed a sample of 20 luminous SDSS 
quasars at z ~ 6, and find that the mean Lya EW and 
FWHM are 56 A and 25 A (also corrected for the Lya 
forest absorption) with large scatters of 40 A and 11 
A, respectively. In Table 2, three of our quasars (SDSS 
J0005-0006, SDSS J0303-0019, and SDSS J2315-0023) 
have Lya FWHM less than half of the typical value. 
The Lya FWHM in both SDSS J0005-0006 and SDSS 
J0303-0019 is only - 1600 km s" 1 . However, their EW 
are close to or stronger than the typical Lya EW. On 
the contrary, the other two quasars, SDSS J0203+0012 
and SDSS J205406.49-000514.8 (hereafter SDSS J2054- 
0005), have typical Lya FWHM, but very weak Lya EW. 

The best-fitting power-law continuum is also used to 
calculate toi 45 o and Mi 45 o, the apparent and absolute 
AB magnitudes of the continuum at rest-frame 1450 A. 
The results are given in Table 3. The sample spans a 
luminosity range of —26.5 < M1450 < —25.4. Because 
the Lya emission usually consists of a large fraction of 
the total emission in the z-band spectra of these quasars, 
the large scatter in the Lya EW results in a large scatter 
in the distributions of mi4so and M1450, even though the 
zab magnitudes lie in the small range 20.5 < zab < 20.9. 

3.1. Notes on individual objects 

SDSS J0005-0006 (z = 5.850). SDSS J0005-0006 
was discovered by Fan et al. (2004). This quasar has a 
very narrow Lya emission line. The rest-frame FWHM 
of Lya is only 1680 km s _1 (after corrected for the Lya 
forest absorption). The strong Nv emission line is well 
separated from Lya. The central BH mass is 3 x 10 8 M Q 
(Kurk ct al. 2007), an order of magnitude lower than the 
BH masses in luminous quasars at z ~ 6 (e.g. Barth et 
al. 2003; Vestergaard 2004; Jiang et al. 2007; Kurk et al. 
2007). SDSS J0005-0006 was marginally detected in the 
Spitzer IRAC 8.0/xm band and was not detected in the 
Spitzer MIPS 24^im band, indicating that there is no 
hot dust emission in this quasar (Jiang et al. 2006a). 

SDSS J0203+0012 (z = 5.854). SDSS J0203+0012 
was independently discovered by matching the UKIDSS 
data to the SDSS data (Venemans et al. 2007). Its Lya 
emission line is broad but weak. The Ly/3 and O vi A1033 
(hereafter O vi) emission lines are clearly seen at ~ 7000 
A. A Civ AA1548,1550 absorption doublet is detected at 
AA = 8480.1, 8494.3 A; this was also noticed by Vene- 
mans et al. (2007). 

SDSS J0303-0019 (z = 6.070). SDSS J0303-0019 
also has a very narrow Lya emission line. The rest-frame 
FWHM of Lya is 1580 km s^ 1 , similar to the Lya width 
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Fig. 4. — Optical spectra of the six high-redshift quasars discovered in the SDSS deep stripe. The spectra of SDSS J0005— 0006, SDSS 
J0353+0104, and SDSS J2315— 0023 were taken on Keck/ESI with a total exposure time of 60 minutes on each source. The spectra of the 
other three quasars were taken on Magellan/LDSS-3 with a total exposure time of 100 minutes on each source. The ESI spectra have been 
binned by 10 pixels and the LDSS-3 spectra been smoothed by 5 pixels. Four quasars in this sample are at z > 6. Three quasars have 
narrow Lya emission lines. SDSS J0353+0104 is a BAL quasar. 



TABLE 2 

Properties of the Lyo and Nv Emission Lines 



Quasar (SDSS) 


Rcdshift 


EW (Lya) 


FWHM (Lya) 




EW (N v) 


FWHM (Nv) 


J0005-0006 


5.850 


81.5±2.5 


6.8±0.4 (1680 km s 




25.0±1.0 


18.9±0.6 


J0203+0012 


5.854 


35.9±1.5 


31.0±2.4 (7650 km 


3-1) 


8.1±0.5 


20.6±0.9 


J0303-0019 


6.070 


139.4±2.4 


6.4±0.5 (1580 km s 




24.4±0.6 


9.9±0.2 


J0353+0104 a 


6.049 








J2054-0005 


6.062 


17.0±1.1 


19.4±3.8 (4890 km 


s-i) 


12.8±0.7 


30.8±1.4 


J2315-0023 


6.117 


126.8±3.2 


9.8±0.5 (2420 km s 


- 1 ) 


37.4±1.4 


17.7±0.7 



Note. — Rest-frame FWHM and EW arc in units of A. The EW and FWHM of the Lya 
emission lines have been corrected for Lya forest absorption. 
a This is a BAL quasar, so we did not measure its emission line properties. 
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TABLE 3 

Continuum Properties of the Quasars 



Quasar (SDSS) 


Redshift 


mi450 (mag) 


M1450 (mag) 


J0005-0006 


5.850 


20.83±0.10 


-25.82±0.10 


J0203+0012 


5.854 


20.94±0.10 


-25.72±0.10 


J0303-0019 


6.070 


21.28±0.07 


-25.43±0.07 


J0353+0104 


6.049 


20.22±0.08 


-26.49±0.08 


J2054-0005 


6.062 


20.60±0.09 


-26.11±0.09 


J2315-0023 


6.117 


21.34±0.08 


-25.38±0.08 



-27 



-26 



r 



-25 



of SDSS J0005-0006. The strong N v emission line is well 
separated from Lya. The Ly/3 and O vi emission lines are 
clearly detected at ~ 7300 A. The continuum emission 
in this quasar is very weak. The absolute magnitude 
^1450 is —25.43, roughly two magnitudes fainter than 
the luminous SDSS quasars at z ~ 6. 

SDSS J0353+0104 (z = 6.049). SDSS J0353+0104 
is a BAL quasar, as seen from strong absorption features 
around the Lya emission line. Its redshift was measured 
from O I. The fraction of BAL quasars in this small 
sample is one out of six, similar to the low-redshift frac- 
tion (Trump et al. 2006), although C IV observations may 
yield further BAL examples. The Ly/3 and O vi emission 
lines are seen at ~ 7300 A. 

SDSS J2054 0005 (z = 6.062). SDSS J2054-0005 
has a very weak Lya emission line. The rest-frame EW of 
Lya is only 17.0 A, significantly smaller than the typical 
EW. But the FWHM of Lya, ~ 30 A, is similar to the 
mean value of Lya FWHM. 

SDSS J2315-0023 (z = 6.117). SDSS J2315-0023 is 
the most distant quasar in this sample. The properties 
of the Lya and N v emission lines are similar to those of 
SDSS J0005-0006 and SDSS J0303-0019. It has a narrow 
but strong Lya emission line. The rest-frame EW and 
FWHM of Lya are 127 A and 2420 km s _1 , respectively. 
It also has a very strong N v emission line. 

4. QLF AT Z ~ 6 

The six quasars presented in this paper provide a flux- 
limited quasar sample at z > 5.8. The survey area is 
260 deg 2 and the magnitude limit is zab =21. In this 
section we calculate the spatial density of the z > 5.8 
quasars in the SDSS deep stripe, and combine this faint 
quasar sample with the SDSS bright quasar sample to 
derive the QLF at z ~ 6. 

We use the selection function to correct the sample 
incompleteness due to the selection criteria we applied. 
The selection function is defined as the probability that 
a quasar with a given magnitude, redshift, and intrin- 
sic spectral energy distribution (SED) meets our selec- 
tion criteria. By assuming a distribution for the intrinsic 
SEDs, we calculate the average selection probability as 
a function of magnitude and redshift. To do this, wc 
first calculate the synthetic distribution of quasar colors 
for a given (Mi 4 5o,z), following the procedures in Fan 
(1999) and Fan et al. (2001a). Then we calculate the 
SDSS magnitudes from the model spectra and incorpo- 
rate photometric errors into each band. For an object 
with given (M 1450 ,z), we generate a database of model 
quasars with the same (M1450, z). The detection proba- 
bility for this quasar is then the fraction of model quasars 
that meet the selection criteria. The details of the model 
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Fig. 5. — Quasar selection function as a function of M1450 and 
2 for the two selection criterion equations (la and lb) based on J 
and H bands. The contours in the figure are selection probabilities 
from 0.9 to 0.1 with an interval of 0.1. The solid circles are the 
locations of the six 2 ~ 6 quasars in our sample. 

and simulation are described in Fan (1999) and Fan et 
al. (2001a). 

Figure 5 shows the selection function as a function of 
M1450 and z for the two selection criteria (Equations 1 
and 2) based on J and H bands. The contours in the 
figure are selection probabilities from 0.9 to 0.1 with an 
interval of 0.1. The sharp decrease of the probability 
at z ~ 5.8 is due to the color cut of iab — zab > 2.2. 
The two selection functions are slightly different. Due 
to smaller photometric errors in the i and z bands, our 
survey probes ~ 1.5 magnitude deeper than the SDSS 
main survey (see Fan et al. 2001a). The solid circles are 
the locations of the six z <~ 6 quasars. All the quasars 
have the selection probabilities greater than 0.5. 

We derive the spatial density of the z > 5.8 quasars us- 
ing the traditional \jV a method (Avni & Bahcall 1980). 
The available volume for a quasar with absolute magni- 
tude M1450 and redshift z in a magnitude bin AM and 
a redshift bin Az is 



V a = ( I p(M 1450 ,z)^-dzdM, 

JAM JAz dz 



(3) 



where p{M\^, z) is the selection function used to cor- 
rect the sample incompleteness. We use one M1450-Z bin 
for our small sample. The redshift integral is over the 
redshift range 5.7 < z < 6.6 and the magnitude integral 
is over the range that the sample covers. The spatial 
density and its statistical uncertainty can be written as 



E^77' a ^= E(^7 



1/2 



(4) 
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Fig. 6. — QLF at z ~ 6. The filled circle represents the density 
of the six quasars in the SDSS deep stripe and the filled triangles 
represent the densities from a study of 17 quasars from the SDSS 
main survey. The open circle is the constraint from the quasar 
discovered by Cool et al. (2006). The dashed line shows the best 
power-law fit to the SDSS quasars. The slope of the QLF is — 3.1 ± 
0.4. 

where the sum is over all quasars in the sample. This 
is similar to the revised 1/V a method of Page & Car- 
rera (2000), since p(M\^Q, z) has already corrected the 
incompleteness at the flux limit. We find that the 
spatial density at (z) — 6.0 and (Mi 45 o) = —25.8 is 
p = (5.0 ± 2.1) x 10~ 6 Mpc- 3 mag" 1 . 

In the SDSS main survey, more than 20 quasars at 
z > 5.8 have been discovered. Seventeen of them were se- 
lected using similar criteria, and consist of a flux-limited 
quasar sample with zab < 20. This bright quasar sample 
includes 14 published quasars (e.g. Fan et al. 2006a) and 
three new quasars (Fan et al. in preparation). We com- 
bine this sample with our sample to derive the QLF at 
z ~ 6. The quasars in the combined sample are divided 
into four luminosity bins as shown in Figure 6. The QLF 
at the bright end at z ~ 6 is well fit to a single power 
law <&(Li45o) cx h\ 



-0 



$(Mi45u) = $*10 



-0.4(/3+l)(M 14 so+26) 



(5) 



where we only consider luminosity dependence and ne- 
glect redshift evolution over our narrow redshift range. 
The best fits are $* = (5.2 ± 1.9) x 10~ 9 Mpc~ 3 mag" 1 
and = -3.1 ± 0.4. The slope = -3.1 ± 0.4 is consis- 
tent with the slope — 3.2 ± 0.7 derived from the luminous 
sample alone by Fan et al. (2004) and with the slope 
> —4.63 (3<r) constrained by the lack of lenses in four 
high-redshift quasars (Richards et al. 2004). 

5. DISCUSSION 

We have derived the QLF at z ~ 6 and found a steep 
slope (/? = -3.1) at the bright end. From the SDSS Data 
Release Three, Richards et al. (2006) showed strong evi- 
dence that the bright-end slope of the QLF significantly 
flattens from z ~ 2 to 4. They found that the slope 
at z ~ 2.0 is —3.1 ± 0.1 and the slope at z ~ 4.1 flat- 
tens to —2.1 ± 0.1. At z ~ 6 the slope steepens again to 
— 3.1±0.4 at a significance level of ~ 2.5c. The flattening 
of the slope at z ~ 4 has been questioned by Fontanot et 
al. (2007), who argued that the spectral index Richards 



et al. (2006) used to correct for sample incompleteness 
was too blue. Hopkins et al. (2007) have analyzed the 
bolometric QLF from multiple surveys and stress that 
the flattening is seen at high significance. The lack of 
flattening claimed by Fontanot et al. (2007) would be 
real only if the distribution of quasar SEDs was redshift- 
dependent, contrary to what is found in most observa- 
tions. Hence, the slope change from z ~ 4 to z ~ 6 
is highly likely to be physical. The steepening of the 
slope at z ~ 6 has important consequences in under- 
standing early BH growth in quasars. Quasar evolution 
at z ~ 6 is limited by the number of e-folding times avail- 
able for BH accretion, therefore the shape of the QLF at 
z ~ 6 puts strong constraints on models of BH growth 
(e.g. Wyithe & Loeb 2003; Hopkins et al. 2005; Volon- 
teri & Rees 2006; Wyithe & Padmanabhan 2006; Li et 
al. 2007), and helps determine whether standard models 
of radiatively efficient Eddington accretion from stellar 
seeds are still allowed, or alternative models of BH birth 
(e.g. from intermediate-mass BHs) and BH accretion 
(super-Eddington or radiatively inefficient) are required 
(e.g. Volonteri & Rees 2006). 

The steepening of the QLF slope also has a strong im- 
pact on the quasar contribution to the ionizing back- 
ground at z ~ 6. The reionization of the universe occurs 
at z = 11 ±4 (Spergel et al. 2007) and ends at z ~ 6 (e.g. 
Fan et al. 2006b). Studies have shown that quasars/AGN 
alone are not likely to ionize the IGM at z ~ 6 (e.g. 
Dijkstra et al. 2004; Meiksin 2005; Willott et al. 2005; 
Douglas et al. 2007; Salvaterra et al. 2007; Srbinovsky & 
Wyithe 2007; Shankar & Mathur 2007). Galaxies proba- 
bly can provide enough photons for the reionization (e.g. 
Yan & Windhorst 2004; Bouwens et al. 2006; Kashikawa 
et al. 2006; McQuinn et al. 2007), however, the individual 
contributions of galaxies and quasars to the reionization 
are not well determined. The galaxy contribution is un- 
certain due to our lack of knowledge of factors such as 
the star-formation rate, the faint-end slope of the galaxy 
luminosity function, and the escape fraction of ionizing 
photons from galaxies (e.g. Bunker et al. 2004; Bouwens 
et al. 2006); while the quasar contribution is poorly con- 
strained due to the lack of the knowledge of the faint end 
of the QLF at z ~ 6. 

We estimate the rate at which quasars emit ionizing 
photons at z ~ 6 from the QLF derived in § 4. Following 
Fan et al. (2001a), we calculate the photon emissivity of 
quasars per unit comoving volume at z ~ 6 as 



No 



q p 
e 1450 n 1450> 



(6) 



where e? 450 



is the quasar emissivity at 1450 A in units of 
erg s -1 Hz" 1 Mpc -3 and nf 450 is the number of ionizing 
photons for a source with a luminosity of 1 erg s" 1 Hz -1 
at 1450 A. We estimate ef 450 from 



fc 1450 



$(M; 



1450^1450 



,dM- 



1450- 



(7) 



where $(Mi4 50 ) is the QLF at z ~ 6 and the integral 
is over the range —30 < Mi 450 < —16. At low redshift, 
the shape of the QLF can be well modeled as a double 
power law with a steep bright end and a flat faint end 
(e.g. Boyle et al. 2000; Richards et al. 2005; Jiang et al. 
2006b). At z ~ 6, we assume a double power-law form 
for the QLF as well, 

$(Mi45o) = 



10 



10 0.4(a+l)(Mi450-M* 460 ) + 10 0.4(/3+l)(M 14 5O-M 1 * 450 ) 



,(8) 



where (3 is the bright-end slope of the QLF, a is the faint- 
end slope, and Af 4450 corresponds to the characteristic 
luminosity L^ 450 . At z < 2, a is roughly —1.6 and Mj* 450 
is a function of rcdshift. Little is known about a and 
M* 450 at z > 4. At z ~ 6, the lower limit of Mj* 450 is 
about —25 as seen from Figure 6. Therefore we calculate 
e i450 f° r a ran g e °f a between —1.2 and —2.2 and for a 
range of Mj* 450 between —21 and —25. The bright-end 
slope (3 is fixed to —3.1 at the first step. We assume 
$* to be the value derived from Equation 5. This is a 
good approximation in the a and M 4450 ranges that we 
investigate here. 

We calculate n 4450 by assuming the quasar SED fol- 
lowing 

/ v-°- 5 , if A > 1050 A; ... 
K 1 „-!■», if A < 1050 A, < 9 > 



and integrating the SED over an energy range of 1-4 
Rydberg. The SED slope at A > 1050 A is taken from 
Vanden Berk et al. (2001) and at A < 1050 A from Zheng 
et al. (1997). 

The total photon emissivity per unit comoving volume 
required to ionize the universe is estimated to be 



Af ion {z) = l0 512 M P c- i s 



l + z 



n b h 2 

0.02 



(10) 

(Madau et al. 1999), where baryon density Qfc/r = 0.02 
(Spergel et al. 2007) and C is the clumping factor of the 
IGM. We examine the effects of three values 1, 10, and 
30 for C. 

The results are shown in Figure 7(a). The solid lines 
represent the photon emissivity per unit comoving vol- 
ume from quasars as a function of the faint-end slope 
a and the characteristic luminosity M 4450 . The dashed 
lines show the photon emissivity required to ionize the 
IGM at z ~ 6 for C = 1, C = 10, and C = 30, respec- 
tively. There is an uncertainty of 0.4 (la) in the bright- 
end slope [3, so we repeat our analysis for (3 = —3.5. The 
results are given in Figure 7(b). This figure shows that 



1 



The significance of the quasar contribution to the 
ionizing background at z ~ 6 depends on a, (3, 



Mf 450 , and C. 



2. For a given (3 (e.g. (3 = —3.1), the photon emissiv- 
ity of quasars is a strong function of -M 4450 , but is 
not sensitive to a as long as a > —2.2. Note that 
a is between —1.0 and —2.0 at low rcdshift. 

3. In Equation 10, Mi on (z) is a linear function of the 
clumping factor C, which is critical here to deter- 
mine whether or not quasars alone can ionize the 
universe. If C = 30 (Gncdin & Ostriker 1997), 
quasars can barely provide enough ionizing pho- 
tons even if M 1 * 450 > —22. For a homogeneous 
IGM (C = 1), however, quasars can easily provide 
enough ionizing photons. 

4. The photon emissivity from quasars increases sig- 
nificantly as (3 steepens from —3.1 to —3.5. 
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Fig. 7. — Photon emissivity (per unit comoving volume) from 
quasars as a function of the faint-end QLF slope a and the char- 
acteristic luminosity Afj" 45 „. The solid lines represent the total 
photon emissivity from quasars. The dashed lines show the pho- 
ton emissivity required to ionize the IGM at z ~ 6 for the clumping 
factors C = 1, C = 10, and C = 30, respectively. The upper panel 
is for the bright-end QLF slope (i = —3.1 and the lower panel is 
for (3 = -3.5. 

Since the total photon emissivity of quasars is not sen- 
sitive to a, the quasar contribution to the reionization 
mainly depends on Mj* 450 and C for a given (3. It is clear 
that the quasar/AGN population can provide enough 
photons required to ionize the IGM only if the IGM is 
very homogeneous or the break luminosity is low. At low 
redshift, Mj 450 varies with redshift and the typical value 
is between -23 and -21. At z ~ 6, the lower limit of 
M 4450 is about —25. To measure M 4450 at z ~ 6, much 
deeper surveys are needed. The clumping factor C is 
also poorly constrained. Previous studies used a range 
of values from 1 to ~100. If the typical value of C is 
10-30 (e.g. Gnedin & Ostriker 1997; Madau et al. 1999), 
from Figure 7(a) the quasar population with (3 = —3.1 
is not likely to provide enough ionizing photons. How- 
ever, if the reionization occurs outside-in and denser gas 
is ionized at a later time when most of the volume of 
the universe has been reionized, dumpiness does not sig- 
nificantly increase the number of photons required for 
reionization (Miralda-Escude et al. 2000). In this case, 
the equivalent C is close to 1, and quasars can provide 
the required number of photons to ionize the IGM. 

6. SUMMARY 

We have discovered five quasars at z > 5.8 in 260 deg 2 
of the SDSS deep stripe, including one previously dis- 
covered by Venemans et al. (2007). The most distant 
one is at z = 6.12. These quasars were selected as i- 
dropout objects from the coadds of 10 SDSS imaging 
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runs, going ~ 1.5 magnitudes fainter than the SDSS 
main survey. The five quasars, with 20 < zab < 21, 
are 1-2 magnitudes fainter than the luminous z ~ 6 
quasars found in the SDSS main survey. The Lya emis- 
sion lines in two quasars SDSS J0303-0019 and SDSS 
J2315-0023 are narrow (FWHM ~ 1600 and 2400 km 
s" 1 ) but strong (EW - 139 and 127 A), while the Lya 
emission lines in another two quasars SDSS J0303-0019 
and SDSS J23 15-0023 are broad (FWHM - 7700 and 
4900 km s" 1 ) but weak (EW ~ 36 and 17 A). The fifth 
one, SDSS J0353+0104, is a BAL quasar. 

The new quasars, together with a previously discovered 
quasar, SDSS J0005-0006, comprise a flux-limited quasar 
sample with zab < 21 at z <~ 6 over 260 deg 2 . The sam- 
ple covers the luminosity range —26.5 < M1450 < —25.4. 
The spatial density of the quasars at (z) = 6.0 and 
(M1450) = -25.8 is (5.0±2.1) x 10~ 9 Mpc" 3 mag" 1 . We 
use a single power-law form to model the bright-end QLF 
at z ~ 6 and find a slope of —3.1 ± 0.4, which is signifi- 
cantly steeper than the slope of the QLF at z ~ 4. Using 
the derived QLF, we find that the quasar/AGN popula- 
tion can provide enough photons required to ionize the 
IGM at z ~ 6 only if the IGM is very homogeneous and 
the characteristic luminosity of the QLF is very low. To 
put better constraints on the quasar contribution, much 
deeper surveys are needed. 

The quasars in this paper were selected from the SDSS 
coadded images with 5-18 runs. Currently the SDSS 
deep stripe has been scanned between 40 and 50 times, 
reaching 2^3 magnitudes deeper than the main survey 
when co-added. We are performing a deeper survey of 
z ~ 6 quasars down to zab ~ 22 in this region. We 
expect to obtain a flux-limited sample with zab < 22 in 
the next few years. 



We acknowledge support from NSF grant AST- 
0307384, a Sloan Research Fellowship and a Packard Fel- 
lowship for Science and Engineering (LJ, XF). We thank 
the MMT staff, Magellan staff, and Keck staff for their 
expert help in preparing and carrying out the observa- 
tions. 

Funding for the SDSS and SDSS-II has been pro- 
vided by the Alfred P. Sloan Foundation, the Partic- 
ipating Institutions, the National Science Foundation, 
the U.S. Department of Energy, the National Aeronau- 
tics and Space Administration, the Japanese Monbuka- 
gakusho, the Max Planck Society, and the Higher Ed- 
ucation Funding Council for England. The SDSS Web 
Site is http://www.sdss.org/. The SDSS is managed by 
the Astrophysical Research Consortium for the Partic- 
ipating Institutions. The Participating Institutions are 
the American Museum of Natural History, Astrophysical 
Institute Potsdam, University of Basel, Cambridge Uni- 
versity, Case Western Reserve University, University of 
Chicago, Drexel University, Fermilab, the Institute for 
Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear As- 
trophysics, the Kavli Institute for Particle Astrophysics 
and Cosmology, the Korean Scientist Group, the Chi- 
nese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astron- 
omy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, Ohio State Univer- 
sity, University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observa- 
tory, and the University of Washington. 



REFERENCES 



Adelman-McCarthy, J. K., et al. 2007a, ApJS, in press (astro- 
ph/0707.3380) 

Adelman-McCarthy, J. K., et al. 2007b, ApJS, submitted (astro- 

ph/0707.3413) 
Avni, Y., & Bahcall, J. N. 1980, ApJ, 235, 694 
Barth, A. J., Martini, P., Nelson, C. H., & Ho, L. C. 2003, ApJ, 

594, L95 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 
Becker, R. H., et al. 2001, AJ, 122, 2850 

Bertin, E., Mellier, Y., Radovich, M., Missonnicr, G., Didelon, P., 
& Morin, B. 2002, Astronomical Data Analysis Software and 
Systems XI, 281, 228 

Bouwcns, R. J., Illingworth, G. D., Blakcslce, J. P., & Franx, M. 

2006, ApJ, 653, 53 

Boyle, B. J., Shanks, T., Croom, S. M., Smith, R. J., Miller, L., 

Loaring, N., & Heymans, C. 2000, MNRAS, 317, 1014 
Bunker, A. J., Stanway, E. R., Ellis, R. S., & McMahon, R. G. 

2004, MNRAS, 355, 374 
Chiu, K., Fan, X., Leggett, S. K., Golimowski, D. A., Zheng, W., 

Geballe, T. R., Schneider, D. P., & Brinkmann, J. 2006, AJ, 131, 

2722 

Cool, R. J., et al. 2006, AJ, 132, 823 

Dijkstra, M., Haiman, Z., & Loeb, A. 2004, ApJ, 613, 646 
Djorgovski, S. G., Castro, S., Stern, D., & Mahabal, A. A. 2001, 
ApJ, 560, L5 

Douglas, L. S., Bremer, M. N., Stanway, E. R., & Lehnert, M. D. 

2007, MNRAS, 376, 1393 

Elston, R. J., et al. 2006, ApJ, 639, 816 
Fan, X. 1999, AJ, 117, 2528 
Fan, X., et al. 2000, AJ, 120, 1167 
Fan, X., et al. 2001a, AJ, 122, 2833 
Fan, X., et al. 2001b, AJ, 121, 31 

Fan, X., Narayanan, V. K., Strauss, M. A., White, R. L., Becker, 

R. H., Pentericci, L., & Rix, H.-W. 2002, AJ, 123, 1247 
Fan, X., et al. 2003, AJ, 125, 1649 
Fan, X., et al. 2004, AJ, 128, 515 
Fan, X., et al. 2006a, AJ, 131, 1203 
Fan, X., et al. 2006b, AJ, 132, 117 



Fan, X., Carilli, C. L., & Keating, B. 2006c, ARA&A, 44, 415 
Fontanot, F., et al. 2007, A&A, 461, 39 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., 

& Schneider, D. P. 1996, AJ, 111,1748 
Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 
Golimowski, D. A., et al. 2004, AJ, 127, 3516 
Goto, T. 2006, MNRAS, 371, 769 
Gunn, J. E., & Peterson, B. A. 1965, ApJ, 142, 1633 
Gunn, J. E., et al. 1998, AJ, 116, 3040 
Gunn, J. E., et al. 2006, AJ, 131, 2332 

Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, J. E. 

2001, AJ, 122, 2129 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Martini, 

P., Robertson, B., & Springel, V. 2005, ApJ, 630, 705 
Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 

731 

Ivezic, Z., et al. 2004, AN, 325, 583 

Jannuzi, B. T., & Dey, A. 1999, in ASP Conf. Ser. 191, Photometric 
Redshifts and High-Rcdshift Galaxies, ed. R. J. Weymann, L. J. 
Storrie-Lombardi, M. Sawicki, & R. J. Brunncr (San Francisco: 
ASP), 111 

Jiang, L., ct al. 2006a, AJ, 132, 2127 

Jiang, L., ct al. 2006b, AJ, 131, 2788 

Jiang, L., ct al. 2007, AJ, 134, 1150 

Kashikawa, N., et al. 2006, ApJ, 648, 7 

Knapp, G. R., et al. 2004, AJ, 127, 3553 

Kurk, J. D. et al. 2007, ApJ, in press (astro-ph/0707.1662) 

Li, Y., et al. 2007, ApJ, 665, 187 

Lupton, R. H., Gunn, J. E., & Szalay, A. S. 1999, AJ, 118, 1406 
Lupton, R. H., Gunn, J. E., Ivezic, Z., Knapp, G. R., Kent, S., & 

Yasuda, N. 2001, in Astronomical Data Analysis Software and 

Systems X, edited by F. R. Harnden Jr., F. A. Primini, and H. 

E. Payne, ASP Conference Proceedings, 238, 269 
Madau, P., Haardt, F., & Rees, M. J. 1999, ApJ, 514, 648 
Mahabal, A., Stern, D., Bogosavljevic, M., Djorgovski, S. G., & 

Thompson, D. 2005, ApJ, 634, L9 
Maiolino, R., Juarez, Y., Mujica, R., Nagar, N. M., & Oliva, E. 

2003, ApJ, 596, L155 



12 



McGrccr, I. D., Becker, R. H., Helfand, D. J., & White, R. L. 2006, 
ApJ, 652, 157 

McQuinn, M., Hernquist, L., Zaldarriaga, M., & Dutta, S. 2007, 

MNRAS, in press (astro-ph/0704.2239) 
Mciksin, A. 2005, MNRAS, 356, 596 

Miralda-Escude, J., Haehnelt, M., & Rees, M. J. 2000, ApJ, 530, 1 

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Page, M. J., & Carrera, F. J. 2000, MNRAS, 311, 433 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, 

S. M., Lupton, R. H., & Ivezic, Z. 2003, AJ, 125, 1559 
Richards, G. T., et al. 2004, AJ, 127, 1305 
Richards, G. T., et al. 2005, MNRAS, 360, 839 
Richards, G. T., et al. 2006, AJ, 131, 2766 

Salvaterra, R., Haardt, F., & Volonteri, M. 2007, MNRAS, 374, 
761 

Schneider, D. P., Schmidt, M., & Gunn, J. E. 1991, AJ, 101, 2004 

Shankar, F., & Mathur, S. 2007, ApJ, 660, 1051 

Sheinis, A. I., Bolte, M., Epps, H. W., Kibrick, R. I., Miller, J. S., 

Radovan, M. V., Bigelow, B. C, & Sutin, B. M. 2002, PASP, 

114, 851 

Shen, Y., et al. 2007, AJ, 133, 2222 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 
Smith, J. A., et al. 2002, AJ, 123, 2121 



Spergel, D. N., et al. 2007, ApJS, 170, 377 

Srbinovsky, J. A., & Wyithe, J. S. B. 2007, MNRAS, 374, 627 

Stoughton, C, et al. 2002, AJ, 123, 485 

Trump, J. R., et al. 2006, ApJS, 165, 1 

Tucker, D., et al. 2006, AN, 327, 821 

Vanden Berk, D. E., et al. 2001, AJ, 122, 549 

Vcncmans, B. P., McMahon, R. G., Warren, S. J., Gonzalcz-Solarcs, 
E. A., Hewett, P. C, Mortlock, D. J., Dye, S., & Sharp, R. G. 
2007, MNRAS, 376, L76 

Vcstcrgaard, M. 2004, ApJ, 601, 676 

Volonteri, M., & Rees, M. J. 2006, ApJ, 650, 669 

Warren, S. J., et al. 2007, MNRAS, 375, 213 

Willott, C. J., Delfosse, X., Forveille, T., Delorme, P., & Gwyn, 

S. D. J. 2005, ApJ, 633, 630 
Willott, C. J., AJ, submitted (astro-ph/0706.0914) 
Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 
Wyithe, J. S. B., & Padmanabhan, T. 2006, MNRAS, 372, 1681 
Yan, H., & Windhorst, R. A. 2004, ApJ, 612, L93 
York, D. G., et al. 2000, AJ, 120, 1579 

Zheng, W., Kriss, G. A., Telfer, R. C, Grimes, J. P., & Davidson, 
A. F. 1997, ApJ, 475, 469 



